Possible Shape Coexistence and Magnetic Dipole Transitions in ^^C and ^^Ne 
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Magnetic dipole(Ml) transitions of N=ll nuclei, ^^C and ^^Ne are investigated by using shell 
model and deformed Skyrme Hartree-Fock+blocked BCS wave functions. Shell model calculations 
predict well observed energy spectra and magnetic dipole transitions in ^^Ne, while the results are 
rather poor to predict these observables in ^'^C. In the deformed HF calculations, the ground states 
of two nuclei are shown to have large prolate deformations close to /32=0.4. It is also pointed out 
that the first K'' = 1/2+ state in ^^Ne is prolately deformed, while the first K'' = 1/2+ state in ^'^C 
is predicted to have a large oblate deformation being close to the ground state in energy. We point 
out that experimentally observed large hindrance of Ml transition between I^ — 1/2+ and 3/2+ in 
^^C can be attributed to a shape coexistence near the ground state of ^^C. 
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I. INTRODUCTION 

Recently, many experimental and theoretical efforts 
have been paid to study structure and reaction mech- 
anism in nuclei near drip lines. It has been known that 
electromagnetic observables give useful information to 
study structure of nuclei, not only ground states but also 
excited states. These observables are expected to pin 
down precise information of the configuration and the de- 
formation in nuclei. Advanced experimental instruments 
reveal several unexpected structure of light nuclei with 
the mass number A~(10-20). One of the current issues 
is a large quenching of magnetic dipole (Ml) transition 
between the first excited 1/2+ state and the ground state 
with /■^=3/2+ in ^^C [l| in comparison with the corre- 
sponding transitions in one of N=ll isotones, ^^Ne [2|. 

The deformation manifests itself in observables like E2 
and Ml moments. In rcf. [|, deformed Skyrme HF+BCS 
calculations were performed to study the evolution of de- 
formations in C and Ne isotopes The calculated electric 
quadrupole moments and magnetic moments were suc- 
cessfully compared with empirical data. It was pointed 
out that the shell occupancy gives the crucial effect on the 
evolution of the deformation of isotope chains. This de- 
formation driving mechanism due to the shell occup ancy 
has been noticed as the nuclear Jahn- Teller effect [10| . 
which gives an intuitive understanding for the evolution 
of deformation. A possible shape coexistence is pointed 
out in ^^C because of different deformation driving effects 
between neutrons and protons. Namely, the first excited 
K"^ = 1/2+ state has oblate deformation and almost de- 
generate with the prolately deformed ground state with 
K'^ = 3/2+. On the other hand, there is no sign of shape 



coexistence in ^^Ne since the shell occupancies are almost 
the same between protons and neutrons. In theoretical 
point of view, it is interesting to see how much differ- 
ences and similarities will appear between the results of 
standard shell model calculations and those of the mean 
field theories. To this end, the HF results are compared 
with shell model results to investigate similarities and 
differences between two models in such observables like 
excitation energies and Ml transitions in -'^'^C and ^^Ne. 

In this paper, we extend the previous calculations in 
ref. 3, and particularly focus on recent experimental 
data of Ml transitions in ^^C and ■^^Nc in order to study 
possible shape coexistence near the ground states of ^^C. 
This paper is organized as follows. We study the energy 
levels, magnetic moments and Ml transitions by using 
shell model wave functions in Section 2. The deformed 
HF-f blocked BCS results are shown in Section 3. A sum- 
mary is given in Section 4. 



II. 



SHELL MODEL CALCULATIONS OF ^^C 

AND 2^NE 



In light and medium mass nuclei, the shell model is one 
of the most successful theories to describe nuclear struc- 
ture in both the ground states and the excited states. 
Shell model calculations are performed in (p-sd) model 
space for ^^C and (sd) model space for ^^Ne with three 
effective interactions PSDMK2 [3, SFO d and WBP [|. 
The excitation energies of the first 1'^= 3/2+, 1/2+ and 
5/2+ states are tabulated in Table |T1 The SFO interac- 
tion is identical to PSDMK2 in (sd) model space so that 
two results are the same for ^-'^Ne. The excitation ener- 



2 



TABLE I: Shell model calculations of excitation energies in 
^'^C and ^^Ne. The shell model calculations are performed by 
using effective interactions PSDMK2, SFO and WBP. For sd 
shell configurations, the interaction matrices of PSDMK2 and 
SFO are the same. Experimental data are taken from ref. [Jl 
for ^'^C and from ref. 2] for ^^Ne. All units are in MeV. 





int. 


3 + 
21 


1 + 

21 


5 + 
21 


5 + 
22 


1 + 

22 




MK2 
SFO 
WBP 


0.305 
0.304 
0.0 


0.0 
0.0 
0.295 


0.711 
0.654 
0.032 


1.679 
1.678 
1.998 


5.007 
5.039 
5.034 




cxp 


0.0 


(0.212) 


0.333 








int. 


3 + 
2 1 


1 + 

2 1 


5 + 
21 


5 + 
22 


1 + 

22 


^^Ne 


MK2 (SFO) 
WBP 


0.0 
0.0 


1.930 
2.870 


0.495 
0.249 


3.250 
3.484 


4.688 
5.815 




cxp 


0.0 


2.794 


0.351 


3.735 





TABLE IL Magnetic moments in ^'^C and ^^Ne in unit of 
The shell model calculations are performed by using effective 
interactions PSDMK2, SFO and WBP with the bare g factors. 
The values in the brackets for '^^ in ^^C and |^ and |^ for 
^^Ne are obtained by using the effective g factors for the IV 
channels, Sqs — —0.2gi^ Tz,Sgi = — 0.15r2 and gp — —I.Otz in 
Eq. ([TJ. Experimental data of magnetic moments are taken 
from ref. le] for ^''C and from ref. flT^ for ^^Ne. 





3 + 


1 + 

21 


5 + 
21 


5 + 
2 2 


1 + 

2 2 


MK2 
SFO 
WBP 


-0.710 (-0.686) 
-0.725 (-0.713) 
-0.858 ( -0.819) 


-1.548 
-1.548 
-1.566 


-1.453 
-1.500 
-1.404 


-1.447 
-1.424 
-1.744 


0.280 
0.232 
0.570 


exp 


10.758(38)1 










2iNe 


3 + 


1 + 

21 




5 + 
22 


1 + 

22 


MK2 
WBP 


-0.887 (-0.720) 
-0.824 (-0.674) 


-1.498 
-1.548 


-0.657 (-0.484) 
-0.643 (-0.518) 


-0.403 
-0.692 


0.228 
0.127 


exp 


-0.661797(5) 




10.49(4)1, 
10.70(8)1, |0.88(20)| 







gies of ^^Nc arc well reproduced by all three shell model 
calculations. It is not surprising since the effective force 
is usually fitted to the data of stable nuclei such as ^^Ne. 
The results of WBP show the best agreement with experi- 
mental excitation energies within lOOkeV difference. The 
calculated results arc much worse in the case of ^^C. The 
interactions PSDMK2 and SFO predict the spin-parity 
of the ground state to be /■^=l/2+, while the observed 
spin-parity is /'^=3/2+. The interaction WBP gives a 
state with 7^=3/2+ as the ground state. However, 5/2+ 
is almost degenerate with the ground state contrary to 
the experimental data. 

Magnetic moments and magnetic dipole (Ml) tran- 
sition probabilities B(M1) are given in Tables and 
IllH rcspectivclv. The magnetic operator is defined as 

f^eff = (ff^^ +<55.)s+ (g?-^ + <5g01 + .9p[n x s](i) (1) 

whereSgs and Sgi are the renormalization factors for the 
spin and the orbital g factors, respectively. The last 
term of Eq. ([1} is the tensor component due to the 
core polarization effect. The shell model results of mag- 



TABLE IIP Shell model B(M1) in and ^^Ne in unit of 
with the bare g factors (the effective g factors). Experimental 
data are taken from ref. [3| for ^'^C and ref. @| for ^^Ne. See 
the captions to Table HIl for details. 
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1+ , 3 + 
2 1 2 1 


5-1- . 3 + 
2 1 2 1 




MK2 
SFO 
WBP 


0.084 (0.045) 
0.077 (0.041) 
0.078 (0.043) 


0.070 (0.031) 
0.077 (0.035) 
0.077 (0.034) 




exp 


O.OlOiO.OOl 


0.082 -f0.032/-0.018 




int. 


1+ 3 + 
21 21 


5+ _^ 3 + 
21 21 


^^Ne 


MK2 
WBP 


0.724 (0.607) 
0.451 (0.390) 


0.173 (0.128) 
0.161 (0.109) 




exp 


0.33±0.05 


0.128±0.03 



netic moments are shown in Tabic [11] with the bare g 
factors and the effective g factors for the IV channels, 
5g, = -0.2gl^T, = -0.2 7° V ^, Jg; = -O.ISt^ and 
gp = — I.Otz. For the magnetic moments, the effect of 
Sgs and 6gi cancel each other largely and that of the 
tensor component gp is very small. The net effect of 
the effective operator is less than 5% in ^^C and 20% in 
^^Ne. In comparison with experimental data, the opti- 
mum quenching factor Sgs depends on the model space 
and the effective interaction. For ^^C, small quenching 
factors {Sgs/gi^ - 0.0 for PSDMK2 and SFO, Sgs/gi^ ~ 
— 0.2tz for WBP) give good agreement with the experi- 
mental data. Slightly larger values {Sgs/gl^ ~ —0.25tz 

for PSDMK2 and SFO, Sg^/gi^ 0.2t^ for WBP) give 

reasonable results in the case of ^^Ne. 

In Table IIIIl two empirical Ml transition probabili- 
ties in ^^Ne are reasonably well reproduced by the shell 
model calculations. The best results among the three 
interaction are given by WBP interaction with the ef- 
fective spin g factor Sgs/gl^ = —0.2tz. We can see 
in Tables [H [U IIIIl that the shell model provides good 
agreement not only for the excitation energies but also 
for the magnetic moments and Ml transition probabili- 
ties in ^^Ne. In -'^^C, the Ml transition probability from 
/'^=5/2+ to 3/2+ is reproduced well by the shell model 
with the bare g factors. However, the transition probabil- 
ity from /'^=l/2+ to 3/2+ is very poorly predicted, i.e., 
the empirical data is almost one order magnitude smaller 
than the shell model predictions with the bare g factors. 
The effective g factors adopted in the magnetic moments 
in Table HIl decrease substantially the B(M1) values in 
^^C. However, these effective g factors do not give any 
satisfactory results for the measured two transitions be- 
tween (r = 5/2+ ^ 3/2+) and (/'^ = 1/2+ ^ 3/2+) as 
shown inside of the brackets in Table IIIIl Recently, the 
description of Ml transitions in ^^C has been consider- 
ably improved with the use of a modified SFO Hamilto- 
nian 0. 
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III. DEFORMATIONS AND MAGNETIC 
DIPOLE TRABSITIONS IN ^^C AND ^^NE 



The neutron number dependence of deformations was 
studied along the chain of C and Ne isotopes in ref. Q by 
performing deformed HF+blockcd BCS calculations with 
Skyrme interactions SGII and SIII. In this study, we per- 
form the same deformed HF calculations of two N=ll 
isotones ^^C and ^^Ne with a different Skyrme interac- 
tion SkO'. We found that the results of SkO' are very 
close to those of SGII and SIII. One advantage of SkO' 
is to give an oblate deformed ground state for ^^C with 
the original spin— orbit interaction , while the spin— orbit 
interaction was reduced in SIII and SGII to obtain the 
oblate deformation. In numerical calculations, the ax- 
ial symmetry is assumed for the HF deformed potential. 
The pairing interaction is taken to be a density dependent 
pairing interaction in BCS approximation For numerical 
details about the pairing calculations, see refs. [1, 

Deformed Skyrme IIF+ blocked BCS results are shown 
in Fig. [1] (a) for "C and Fig. [1] (b) for ^iNe. The defor- 
mation and the intrinsic Qo moments are tabulated in 
Table HVl for ^^C, and ^^Ne. The ground states are pre- 
dicted to be K'^ = 3/2+ state in both nuclei having large 
prolate deformations /?2 =0.366 for ^^C and 0.391 for 
^^Ne, respectively. The spin-parity of calculated results 
can be compared with the observed ones I'^ = 3/2+ in 
both nuclei. In ^^C, the first excited state is predicted 
to be K'^ = 1/2+ state with a large oblate deformation 
P2 = —0.270. The energy difference from the ground 
state is rather small with Ex=0.56MeV. On the other 
hand, the first excited = 1/2+ in ^^Nc is predicted 
to have a large prolate deformation =0.287 with a 
large excitation energy Ex=2.33MeV. This difference in 
K"^ = 1/2+ state can be understood as a manifestation 
of the nuclear Jahn- Teller effect due to the proton config- 
uration [13]. In general, a few particles top of the closed 
shell drives prolate deformation, while a few holes prefer 
oblate deformation. There is a strong competition be- 
tween prolate driving N=ll neutrons and oblate driving 
Z=6 protons in "'^^C. Namely, the Z=6 proton configu- 
ration, two proton holes in the Z=8 closed shell, prefers 
the oblate deformation as is the case of the ground state 
of ^^C, while the N=ll neutrons tends to drive prolate 
deformation. Consequently, in ^^C, the ground state is 
prolately deformed due to the effect of neutron config- 
uration. However, the first excited = 1/2+ state 
becomes oblate under the influence of the deformation 
driving force of protons. In ^^Ne, both the proton and 
neutron configurations drive prolate deformation so that 
there is no sign of the shape coexistence. The observed 
excitation energy of the first I'^ = 1/2+ state is very low 
in I'^C as Ex=0.212MeV, while that of ^^Ne is higher as 
Ex=2.79MeV. These experimental observations are con- 
sistent with the calculated results in Tabic HVl as far as the 
excitation energies are concerned. Thus we identify the 
first excited I"" = 1/2+ state as K'' = 1/2+ in both ^^C 
and ^^Ne with different large deformations P2 = —0.270 



and 0.287, respectively. The / = 1/2+ in ^^Ne was inter- 
preted in [III as the head of rotational band with a large 
prolate deformation [ll[. The = 5/2+ state is ob- 
served at very low excitation energy around Ex=0.3MeV 
in both nuclei. In the HF calculations, no K'^ = 5/2+ 
state appears at the energy below Ex ^ IMcV. Thus, we 
interpret that the observed first excited = 5/2+ state 
in both nuclei is a member of the rotational band with 
K"^ = 3/2+. In the case of ^^Ne, the ground state and 
the first excited state were identified as members of the 
same rotational band giving consistent predictions of as- 
sociated observed properties 1^. This interpretation is 



also supported by the large deformation length observed 
in the excitation to I'' = 5/2+ state in the proton inelas- 
tic scattering on -'^^C [l3| . 
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FIG. 1: Energy surfaces as a function of deformation param- 
eter /32 in "C and ^^Ne. Deformed HF-|-blocked BCS calcu- 
lations are performed with a Skyrme interaction SkO'. 

We study the magnetic dipole transitions between the 
excited and ground states in ^''C and ^^Ne using the de- 
formed HF wave functions. 

For axially symmetric deformation, the deformed 
many-particle initial and final states are expressed as 
a direct product of neutron and proton single-particle 
states; 



l^) = k)k), 



(2) 
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TABLE IV: Energies, deformations, Q moments and magnetic 
moments in ^'^C and ^^Ne with a Skyrme interaction SkO'. 
The magnetic moment ^ is calculated for I — K state with 
the bare neutron g— factor. Experimental data are the same 
as for Table HTl (experimental uncertainties are omitted). 



The transition matrix element can be written for one 
neutron quasiparticle states as 





A"" 


Ex 


ft 


Qop 
(fm^) 


Qon 

(fm^) 
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^i(exp) 
(Miv) 


17p 


3 + 
2?, 


u.u 


U.ODD 


1 9/1 


oo.uo 


1 1 Q7 


-U.o / / 


|U. ( Oo| 




0.56 


-0.270 


-15.27 


-35.94 


-3.420 


-1.767 






0.99 


-0.247 


-13.40 


-34.43 


-0.764 


-1.126 






1.21 


0.272 


13.59 


40.75 


-1.101 


-0.947 




2iNe 


3 + 

!^ 

5 + 
21 

1 + 
22 


0.0 


0.391 


42.15 


46.98 


-1.112 


-0.728 


-0.661797 




2.33 


0.287 


31.67 


32.05 


-2.557 


-1.523 






2.92 
4.91 


0.226 
0.357 


25.85 
39.570 


23.83 
43.547 


-0.764 
-0.231 


-1.040 
-0.594 


|0.49| ~ 
|0.88| 



{v'TT'K'\M{Ml)\iynK) = 
{,.'\ax,K'iP'2)MiMl)„aiM)W)W\^) 



(5) 



where (tt'Itt) and {i^'liy) are quasi-particle vacuum over- 
laps of neutrons and protons, respectively. 

The in-band Ml transition probability can be written, 



for a band with K > ^, a,s 



B{Ml-Kh^K,h^h±^) 



3 ,,2 



(6) 



TABLE V: Ml transition probabilities B(M1) in ^''C and 
2^Ne in unit of ^%[. The deformed HF calculations are per- 



where (?_r = ^ and gx is the intrinsic g factor, KgK = 
(A'|g;?3 + ggS^\K) . The magnetic moment is expressed as 

if2 



M = 9rI + (.9 A' - gn) 



1 



(7) 



formed by using a Skyrme interaction SkO'. The effective spin 
g-factor glff = 0.5g*'"''=(0.7g^'''=) is adopted. Experimental 
data are taken from ref. [l| for ^^C and from ref. [2| for ^^Ne. 





/l(i^) = i+-J2(J<') = f: 






0.00068(0.00094) 


0.116(0.179) 


exp 


O.OlOiO.OOl 


0.082 -h0.032/-0.018 whprP 


2^Nc 


0.208(0.272) 


0.152(0.224) 


exp 


0.33±0.05 


0.128±0.003 G 



where the component of the total angular momen- 
tum along the symmetry axis is denoted by K [l3 . 
[Tsl and — a'^p^a^p^- ■ ■\f{l32)) denotes the multi- 

quasiparticle neutron (proton) state. The state |/(/32)) 
is the quasiparticle vacuum with deformation (32- The 
quasiparticle operator is connected to the real particle 
operators c^{P2) and c(/32) in the deformed basis by 



A/j, ^ 



(3) 



where A specifies the quantum numbers of Nilsson orbit, 
v\p.{(32) is the BCS occupation amplitude and u\p{/32) ~ 
Vl - vxp{P2)'^. The operators c^(/32) and CAp(/32) are 
further expanded by the spherical bases as c\^{f32) = 

'^Aa(/^2)cj,^ whcrc the amplitude d^^[j32) is denoted 
by the quantum numbers a = (n, /, j). 

The intrinsic Ml single-particle transition operator is 
expressed as 

M{M\) = X 
(E. - 9r) li + - 9r) Si) + 9rI) (4) 

where g/(i), gs{i) and g^i are the orbital, spin g factors 
and the gyromagnetic ratio of the rotor respectively, in 
unit of the nuclear magneton fiN — eh/2mpC. 



For K ^ K' case, the Ml transition probability is written 
to be 

B{M\;Kh^K'a2)^i^[i%y~ 
<hKlK' -K\l2K' >^ G'' <t:'\t^>^ (8) 



+ (.9/ - 9r) {^'\ayK'{l3'2)lAKa{K{P2)\y)- (9) 

with /^K = K' - K. 

The calculated B(M1) values are tabulated in Table 
El We adopt g^Jf = O.^g^"'''' and g^Jf = O.lg''/'''' for the 
effective neutron spin g-factor in the calculations. One 
can see a large hindrance in B(M1) from Ii = 1/2+ to 
Ii — 3/2+ transition in ^^C. This is entirely due to the 
shape difference between the ground and the first ex- 
cited states. Namely, the value of the core overlap of 
BCS vacuums < q'\q >=< >< tt'Itt > is calculated 
to be < g'(/32 = 0.366)|g(/32 = -0.270) >=0.0378 be- 
tween / ^ 3/2^ and / = l/2j^ states in ^^C. On the 
other hand the corresponding overlap in ^^Ne is close to 
1.0, i.e., < q'iP2 = 0.391)|g(/32 = 0.287) >=0.982 since 
both the initial and the final state have large prolate de- 
formations. For in-band transition (/i = 5/2+ I2 = 
3/2+, K = 3/2), the calculated B(M1) values with the 
g-factor g'^^f ~ O.bgg""^^ agree well with the observed 
ones within the experimental accuracies. This quenching 
factor is somewhat smaller than the adopted values in 
rare-earth nuclei, but it is still in the acceptable range. 
The g-factor g'^^-^ = 0.7 g'^"'^'^ gives better results for the 
transition {h = Ki = 1/2+ ^ h = K2 = 3/2+) in 
^^Ne. The deformed HF calculations gives the lowest 
Ki =1/2+ state having asymptotic quantum numbers 
[iVngAfJ] = [2201/2]. There is another K'' = 1/2+ state 
with [Nn^Ml] = [2111/2] having slightly higher energy. 
The latter has about 2 times larger B{M\) value for 



5 



the transition to the X = 3/2 ground state. It is ex- 
pected that a small mixing of the [NnsAH] = [2111/2] 
state increases the B{M1) value between the (/i ~ Ki = 
1/2+ ^72 = ^2 = 3/2+). We notice that the opti- 
mal quenching spin g factor of the deformed HF results 
is slightly smaller than that of the shell model calcula- 
tions. It is an interesting open question to compare more 
systematically the transition strength of two models in 
a quantitative level. It was mentioned in ref. [l| that 
the halo effect of ^^C may play a role to decrease B(M1) 
value of {I I = 1/2+ I 2 = 3/2+). However no serious 
attempt has been made so far to take into account the 
halo effect on the Ml transitions in ^^C. 

The calculated magnetic moments ^ are shown in Ta- 
ble IIVI The observed magnetic moments show a small 
quenching effect in comparison with the calculated val- 
ues for the ground states of ^^C and ^^Ne. The results 
of deformed HF calculations provide similar quantitative 
predictions to those of the shell models as far as the mag- 
netic moments of the ground states are concerned. The 
observed magnetic moment for the excited I'^ = 5/2+ 
state in ^^Ne is still not accurate enough to perform pre- 
cise comparison with the calculated results. 

The second K = 1/2+ state is found in Table HVl at 
rather low energy =1.21MeV ^^C by the deformed HF 
model, while the / = 1/2^ state is located at '^SMeV 
in the shell model calculations in Table HI So far the 
second 1/2+ is not identified experimentally [l^- It is 
quite interesting to find the 1/2^ state experimentally to 
disentangle the applicability of the two models. 

IV. SUMMARY 

We have studied the magnetic dipole transitions in ^^C 
and ^^Ne using microscopic shell model wave functions 
and deformed HF wave functions. The energy spectra as 



well as Ml transition probabilities of ^^Ne are well repro- 
duced by the shell model calculations, while we need a 
quenching factor for the spin g factor to obtain reasonable 
quantitative agreement. On the other hand, the observed 
Ml transition probability from the first excited 1/2+ to 
the 3/2+ ground state in ^^C was found to be hindered 
by one order of magnitude compared with the shell model 
calculations. The shell model prediction of energy spec- 
tra is also poor in ^''C compared with the experimental 
data. The deformed HF-|-blocked BCS calculations are 
performed with a Skyrme interaction SkO'. The ground 
states of ^^C and ^^Ne are predicted as largely prolate 
deformed states with K'' = 3/2+. In ^iNe, the first 
K"" = 1/2+ state appears at the energy i?x ~2.3MeV 
with a large prolate deformation. On the other hand, 
the first K'^ = 1/2+ state in ^^C has a large oblate 
deformation with E'x '^0.5McV as the result of com- 
petition between the deformation driving force of pro- 
tons and neutrons. The calculated energy difference be- 
tween K'' = 3/2+ and K"" = 1/2+ states is close to 
the observed energy difference between /'^ = 3/2+ and 
r = 1/2+ states both in ^iNe and ^"^0. The strong hin- 
drance of the B(M1) transition between the first excited 
1/2+ to the ground state 3/2+ of ^^C can be attributed 
to the shape difference between the lowest K"" = 1/2+ 
and the first K"" = 3/2+ state as is predicted by the 
deformed HF-|-blocked BCS calculations. 
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